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Abstract 

Aims. We investigate in detail two physical properties of the metric /(if) theory developed by 
Sobouti (2007). We first look for the possibility of producing gravitational waves that travel at the 
speed of light. We then check the possibility of producing extra bending in the lenses produced 
by the theory. 

Methods. We do this by using standard weak field approximations to the gravitational field equa- 
tions that appear in Sobouti's theory. 

Results. We show in this article that the metric theory of gravitation proposed by Sobouti (2007) 
predicts the existence of gravitational waves travelling at the speed of light in vacuum. In fact, 
this is proved in general terms for all metric theories of gravity which can be expressed as powers 
of Ricci's scalar. We also show that an extra additional lensing as compared to the one predicted 
by standard general relativity is produced. 

Conclusions. These two points are generally considered to be of crucial importance in the devel- 
opment of relativistic theories of gravity that could provide an alternative description to the dark 
matter paradigm. 
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1. Introduction 



One of the greatest c hallenges of modern astrophysics is the validation of the dark matter paradigm 
( Bertone et alj|2005 ). Postulating the existence of non-barionic dark matter has given a lot of suc- 
cess to many astrophysical theories. However, no matter h ow hard the strange matter has been 



It was 



looked for , it has never bee n directly observed, nor detected dMunoz]|2004 



Coolevf2006). 



Milg roml (I1983I) who proposed that, to understand certain astrophysical observations 



it was necessary to change Newton's law of gravitation. With time, this idea has developed 
strongly up to the point of building a relativistic Tensor- Vector-Scalar (TeVeS) theory t hat gener- 



alises and substa ntiates the Modified Newtonian Dynamics (MOND) ideas introduced bv lMilgrom 
( Bekensteinll2004 ). 
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The complications introduced by TeVeS have lead different groups to think of an alternative 
possibility. This has been motivated by recent develo pment on metric theories of gravity applied 



to the problem of dark energy. Some researchers (cf . |C apozziello 



Nojiri & Odintsov 



Noiiri & Odintsov 



2003; 



Carroll et al 



2004 



Capozziello & Troisi 



2002; 



2005 



Capozziello et al 



Capozziello et al. 



2003 



2006 



2006, and references therein) have shown that it is possible to explain different 



cosmological observations without the need of dark energy. The idea is to introduce a general 
function f(R) in the Einstein-Hilbert action, instead of the standard Ricci scalar R. The resulting 
differential equations that appear due to this introduction are of the fourth order. This introduces a 
degree of complexity on the field equations, but makes it possible to reproduce some results that 
are usually thought of as being due to a mysterious dark energy field. 



Sobouti 



d2007l) have developed two differ- 



In the same sense, ICapozziello et al.l ((2006) and 
ent f(R) theories that can repr oduce th e anomalous rotation curves produced in different spi- 
ral galaxies. The advantages of ISoboutif s description are many. His theory reproduces naturally 
the standard Tully-Fisher relation, it converges to a version of MOND and, due to the way the 
theory is developed, the resu lting differential equations are of the second order. More impor- 
tantly, |Ca£ozziejlo_etalJ (120061) showed that f(R) = R", with 1. 34 < n < 2.41 reproduces ro- 



tation curves of a number of spiral galaxies. On the other hand, ISoboutil d2007l) showed that if 
f(R) = R^- a / 2 \ then different rotation curves associated to spiral galaxies can be accounted for. 
Also, since »« 1, the modification can be thought of as a small deviation to the Einstein-Hilbert 
action, i.e. f(R) ~ R [1 - (a/2)lnR + (a/2) In (3a)]. 

Central to the development of a good modified theory of gravity that can describe the phe- 
nomenology usually ascribed to dark matter, is the analysis of the propagation of gravitational 
waves and the amount of lensing implied by the theory. In this article, we show that all f(R) = R" 
metric theories of gravity produce gravitational waves that propagate at the velocity of light c in 
vacuum. This is a crucial step in order to consider Sobouti's theory as a possible alternative to the 
dark matter problem and so, it can begin to be applied to real a strophysical situati ons. Relativistic 



Bekensteinl d2006h and references 



theories of MOND have been proposed in the past (see e.g. 
therein) that show superluminal propagation of the waves produced by the fields and so, they were 
rejected immediately. For example, one of the crucial steps in using TeVeS as an alternative to dark 
matter is that the waves produced by the theory are never superluminal. We also show that the the- 
ory proposed by Sobouti describes additional lensing to the standard general relativistic version. 
In other words, the theory developed by Sobouti can in principle be considered an alternative the- 
ory, in order to do astrophysical comparisons with current models of dark matter. Also, this theory 
may challenge what TeVeS has been trying to explain in certain astrophysical situations and so, 
astrophysical predictions between TeVeS and Sobouti's theory must be done in the future. 



2. Modified field equations 

The alternative gravitational model used in w hat follows is o ne that introduces a modification in 
the Einstein-Hilbert action as follows (see e.g. lSobou ti (2007) and references therein): 



f(R)+L m \ ^d 4 x 



(1) 
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where L m is the Lagrangian density of matter and f(R) is an unknown function of the Ricci scalar 
R. Variation of the acti on S with respect to the metric g^ v gives the following field equations 
( Capozzielloetal.12003 ) 



(2) 



where F :- df/dR and is the stress-energy tensor associated to the Lagrangian density of 
matter L m and R^ Y is the Ricci tensor. 

In order to apply to galactic system s, the metric is chosen as a Schwarzschild-like one given 



by dCognola et al 



2005 



Sobouti 2007) 



ds z 



-B(r)dt 2 + A(r)dr 2 + r 2 {dO 2 + sin 2 6 dip 2 ) . 



(3) 



Soboutil showed that the combination A{r)B{r) = g(r) and so g(r) = (r/s) a * 1 + aln(r/s), 



with the Schwarzschild radius s :- 2GM. His calculations also show that the functions A(r) and 
B(r) applied to galactic phenomena are given by 



1 



1 



A (I -a) 
B = 



1 



(;)' 



(l-er/2) 







r\ a 1 
A' 



(4) 
(5) 



For the sake of simp licity we think of f(R) as given by ISoboutiF s (12007h model (but see 



Noiiri & Odintsovi d2004l) for the first model that introduces a logarithmic f(R) in cosmology), 



i.e., 



f(R) = R 



(l-ar/2) 



R 



a a 
1 - -\nR + -ln(3a) 

2 2 



(6) 



The parameter a is chosen in such a way that F(r, a) — > 1, which corresponds to standard general 
relativity, as a — > 0. 



3. Gravitational waves 

Just as it happens in standard general relativity, it is expected that a modified metric theory of 
gravity predicts gravitational waves. These should propagate through space-time with velocity 
equal to that of light. We now show that this happens for all cases in which f(R) = R", where n is 
any number. To do so, we consider a space-time manifold with a metric g^ v deviating by a small 
amount h^ v from the Minkowski metric 77^ in such a way that 

g,iv = V + V' ( ? ) 

with I h^y |<K 1. Using this we can make arbitrary transformations of the coordinates (or refer- 
ence frame) in such a way that jc^ — > xf with ^ small. As usual, we impose the Lorentz gauge 
condition: 



v , M = 0, 



(8) 
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with W v := hf" - (1/2) h rf y and h := h£. The Ricci tensor R MV and the Ricci scalar R to first order 
in h^y are consequently given by 

Rpv * \[^ avJ3 " " v - ^ v ,„ B + i^ v T a <*], (9) 
/?*i[2T' ay , OT + (10) 

With the substitution on equations (0 and (fTOt in the field equations (O with /(/?) = R" we obtain 



6 \ " / (11) 

+ (b - 1)(» - 2)CF,„°)- 2 (T, a ). > ,(T, a ).v. 

Because of the fact that *P = rf*W K „, the gauge condition ([8} and the commutativity of partial 
derivatives to first order in h^ v , the right hand side of equation ([TO is null. Therefore we obtain the 
ordinary wave equation, i.e. *¥p va ' a = 0. 

We now consider a weak perturbation relative to an arbitrary metric ^g^. Then, the metric g^ 
takes the form 

fcnr^pv + V. (12) 

In this expression, | h^ Y | a 0(A/L) where the wavelength /I is small compared to an arbitrary 
characteristic length L, which is related to curvature of space-time. 

Let us define the tensor TV :— - (0) h g^ v . We impose the same transverse traceless gauge 
condition ( Landau & LifshitJl994 ) as we did for equation (0, by replacing the standard derivative 



with a covariant one, i.e. 

W ;y = 0. (13) 

The corrections to the Ricci tensor and the Ricci scalar to first order in h MV are respectively 
given by 

R/3v * ( %v + TV;</ - (a K P , v h a », (14) 

R *® ) R+(% a W la . (15) 

The terms ' ^Rp v , ^R a pfiv an d ^°R are calculated with respect to the metric *• Substitution of 
equations ( fT4l i and (TT3T > in the field equations with the function /(/?) = R n , we find to first order of 
approximation the following relation: 



1 ;/3v 



\vp. a - (0, ^ v » -d( ( V 

+c 2 {(%)- 1 [ % ff ^] ;/3v - (%r 2 [ %«^][ %] . 

-C 3 {(^)" 2 [ ;y [ (0 v^] ;j8 + (%)- 2 [ ,- 



(16) 
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where C\, C2 and C3 are constants for a fixed value of n. The terms involving (0 Rp v , ^Ra/Sfiy and 
(1> R as well as the ones involving first partial derivatives in *P can be neglected due to the fact that 
A/L <K 1 . This result becomes clearer if we propose a solution of equation ( fToT l as 

T /Sv = Re(A /Sve W ) ), (17) 

where the function (p is the eikonal. For the particular case we are dealing with, the eikonal is large 
if we are to satisfy the condition A/L «: 1. 
We now define the 4-wavevector as 

k a := 4>;a. (18) 

Substitution of equation ( fT71 ) in ([Tol l and keeping the dominant terms to second order in </> we find 
the following relation: 

k a k a = 0, (19) 

i.e. the 4-wavevector is null. Therefore, gravitational waves in a f(R) = R" metric theory of gravity 
propagate at the speed of light c. 



4. Gravitational lensing 

A key part of gravitational lensing is the light bending angle due to the gravitational field of a 
point-like mass. This can be determined by using the fact that light rays move along null geodesies, 
i.e. ds 2 = 0. Similarly, its trajectory must satisfy: 

dx^ dx v 

^--=0, (20) 

where £ is a parameter along the light ray. We now derive an expression for the bending angle by 
a static spherically symmetric body. The metric of the corresponding space-time is given by the 
Schwarzschild-like metric (01. Because of the spherical symmetry, the geodesies of d20i > lie in a 
plane, say the equatorial plane 8 = n/2. The fact that the metric coefficients do not depend neither 
on ip nor t, yields the following equations: 

r2 d^ dt = J_ 

df d£ B(r) 

where J is a constant of integration. Substitution of equations ([3]) and (l2Tb into ( f20b and replacing 
£ by ip as an independent variable, we obtain 

1 J 2 (dr\ 2 J 2 



W)-^ A{r) \T v ) -7^°< (22) 

If the closest approach to the lens occurs at a distance r m with an angle <p m , such that (dr m /d<p) = 
0, then J is given by 



(23) 



Equations (l2Tb and ( 122b then yield 
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dip 



(AB) l/2 dr 



(24) 



r^(r/r m ) 2 B(r m )-B( r ) 

We now consider a light ray that originates in the asymptotically flat region of space-time and 
is deflected by a body before arriving at an observer in the flat region. Therefore, equation d24l > 
yields the following expression for the bending angle B: 



B = 2 



(AB) 1/2 dr 



r^(r/r m ) 2 B(r m )-B{r) 



(25) 



To compute the light bending angle in the ISobou ti (2007) f(R) gravitation, we substitute the 



metric functions (O and (0 into expression ((25). If we now define x := r m /r, the deflection angle 
takes the following form: 



B = 2 



f 

Jo 



{sxlr m y (al2) d* 



^B(l)-x 2 B(x) 

This integral can be solved exactly if we assume that sjr m <K 1, in order to obtain: 



(26) 



- .r t ■!(•) r 

JO X"' 2 Vl - X 2 ~ a 2 \ r mJ Jo 



l-a/2 _j (! _ ^3"3«/2j x - a /2 



(1 - x 2 - a ) 



2-a\V 2 



dx } - 7i, 



2 vr 



2-a 



- 1 



+ 2 Vl -a\ — 

' m 



l-a/2 



(27) 



For a = 0, this expression converges to the bending angle expected in traditional general relativity. 
Figure Q] shows the fluctuation AB given by 



A/3 = B-B E 
Be B e ' 



(28) 



where Be is the bending angle obtained by general relativity. This fluctuation is a function of the 
parameter a. In fact, to O(a) it follows that 



&B 1 / r m 

— = -asm 1 

B ¥ . 2 I s 



(29) 



The Figure shows that in order to obtain significant bending, the parameter s/r m needs to be 
not too small for an appropriate value of a. Of course, when s/r m is closer to 1, more bending is 
expected. The result obtained in equation d27j > for the bending angle is only valid for s/r,„ <K 1 and 
so, the plot cannot be used to test greater values of s/r m . However, the trend seen in going from 
s/r m = 10~ 6 to 10~ 2 is strongly suggestive of significant AB/Be for yet smaller impact parameters. 
This could in principle account for anomalous lensing in clusters of gala xies con sidering that the 
releva nt a at those scales might differ from the galactic value calculated by 
(2007) showed empirically that 



Sobouti. In fact. 



Sobouti 



a — ao 



M 



1/2 



(30) 



with ao ~ 2.8 x 10 12 . Even if one assumes that ao is a universal constant, then for very massive 
bodies it may be possible to obtain the required extra bending. 
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0.01 
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Figure 1. The plot shows fluctuations A/3//3e for the bending angle /? as compared to /3e, the one 
predicted by general relativity. From left to right, each plot corresponds to values of s/r m given 
by 10~ 6 , 10~ 5 , 10~ 4 , 1CT 3 , 10~ 2 , respectively. For greater values of s/r m it is expected that more 
bending will be produced. 

T 









- "W 

M 













Figure 2. Diagram of the lensing geometry, i? is the angular position of the source; 9 is the angular 
position of an image; and Dol,Dos and flu are the observe-lens, observer-source and lens-source 
distances, respectively 



5. Lensing framework 

To study how a body acts as a gravitational lens, we begin with a through analysis of lensing by 
a static, spherically symmetric body with mass M. Additionally, we assume the gravitational field 
produced by the lens to be weak. Figure [2] gives a diagram of the lensing situation and defines 
standard quantities: § is the angular position of the source, 8 is the angular position of an image 
and Dol, Dos together with Dis are the observer-to-lens, observer-to-source and lens-to-source 
distances respectively. 



8 S. Mendoza and Y.M. Rosas-Guevara: Waves and lensing in Sobouti's f(R) metric theory 

From the figure, elementary trigonometry establishes the lens equation: 



& = o-p, 



(31) 



where ft :=f3D LS /D s- 

Since r m = 9Dql and using equation OTb . it then follows that the lens equation takes the form 



where 



and 



7T D 



C, := — 



SL 



E D c 



® 2 -®[& + Ci]-C 2 ® a/2 = o 



2 VT 1 ^ 

2-a 



E \2D LS 



-a/2 



(32) 



AGMDls 



® := —, 

E 



:= — . 

E 



(33) 



Dos Doe 

The quantities and <P can be thought of as "scaled angles" with respect to 9e- 

In order to solve the lens equation note that, to first order of approximation, the solution can be 
written as 



= ©o + 01 



(34) 



where 0o represents the standard image position and 0i is the correction term to first order. 
Substitution of equation ( |34l l on ( 1321 gives: 



0, = 



Ci0 o + C 2 ® a J 2 - 1 



2©o -f-d -(a/2)C 2 ® 



a/2-1 



(35) 



We can now obtain the magnification /j :- 6&6I§&& of a lensed image at angular position 6. 
Using equations d34l > and d35l > it follows that this magnification is given by 



_ 0od0o 0i d0 o 



(36) 



With the known positions and magnifications of the images it is now straightforward to obtain 
the time delay At of a light signal. This time is defined as the difference between the light travel 
time r for an actual light ray and the travelled time r eu for an unlensed one, i.e. 



5t — t — T e 



(37) 



To compute the time delay, we use the weak field approximation for the metric, that is goo ~ 
1 + 2(p. For null geodesies it then follows that df = (1 - (p)dl, where / is the Euclidean length. 
Integrating along to the light ray trajectory and introducing the angular variables 9 and § we find 



6t 



DosDol 
D rs 



1 



1 



-{0- §f - -0 E \n(6/e E ) + a0 z HO/Oe) 



(38) 



According to this, if a = 0, we obtain the general relativistic time delay. The extra term 
aO 2 Iii(6/6e) results in a contribution for it when 9 is greater than 9e- For the case of 9 = Oe, 
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this contribution is null. However, the Einstein angle is not a solution of the lens equation d32l even 
if § — 0, so this case is never achieved. In other words, equation (l38b means that the modified time 
delay is greater than what is obtained in standard general relativity. 



6. Conclusion 



We have shown that all metric theories of gravity of the form f(R) = R" produce grav itational 
waves propagating at the velocity of light in vacuum. In part icular, the theory developed by 



Sobouti 



(2007) satisfies this condition. We have also proved that lSoboutir s theory produces an additional 
amount of lensing as compared to standard general relativity calculations. 

Of cour se, more investigation on the physical and astrophysical side of the theory developed 



by ISoboutil needs to be done. By no means can his theory be taken as a fundamental one, but 



rather as a suitable candidate approximation at a certain scale. Particularly, more values of his a 
parameter need to be calculated for different astronomical environments. Also, more development 
in the theory of gravitational lensing needs to be done, in order to compare directly with current 
astronomical data, particularly the anomalou s lensing observed in cosmology. 



As a final remark, we briefly mentio n that 



Sobouti 



s theo ry is not affecte d by the no-go theorem 



proposed by ISoussa & Woodardl d2003l) and expanded by ISoussal (I2003D . This is easy to see if 
we note that some of the components of the perturbations made to the field equations (f2| do not 
necessarily scale with h 2 := W . In fact, in the weak field approximation, the time-component 
of equation (f2]i takes the form 



4 , 

-V 2 d> -2a 

3 w 



In (3a) 



V 2 0+ ^V 2 ln(~2V 2 0) = 



\6nGp 



(39) 



3 V / 3 

From this equation and the value of a given by equation (l30t it is clear that the gravitational 
potenti al 4> may scale as y/GM without ever reaching terms of the order of h 2 . This is the reason 
of why lSoboutir s theory does not satisfy the conditions of the no-go theorem and so, may account 
for an extra amount of lensing as explained before. 

While this article was be ing refereed, and since its appearance in the arXiv, a short c omment 



Sobouti 



was m ade by[Saffari (2007). This author states that some of the calculations made by 
( 12007b are wrong and this is reflected in the circular velocity with a change of sign in one of its 
terms. In his conclusions he states that, although the main results obtained by Sobouti are not 
strongly changed, they may affect higher order corrections of the theory and in the limit applicable 
to very compact objects. Although our article does not reproduce Sobouti's equations mentioned 
by Saffari, we were very well aware of the small error made by Sob outi. In fact, the corre ctions 
to Sobouti's calculations have been published in the Thesis made by iRosas-Guev ara (2006) (see 



equations (4.17), (4.20)-(4.24) of Rosas-Guevara Thesis and compare them with equations (2)-(8) 
in Saffari's article). So, all the results discussed in the present article are not affected in any manner 
by the small error made by Sobouti, because the corrections were already included. 

Although Sobouti's theory seems very attractive it has a small caveat which was discussed by 
him: 'Actions are ordinarily form invariant under the changes in sources. Mass dependence of a (cf. 
equation d30l l) destroys this feature and the claim for the action-based theory should be qualified 
with such reservation in mind". While this fact is not enough to rule out the theory proposed by 
Sobouti, one must not forget it. 
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In summary, it all seems that lSoboutif s theory may play the role of a good candidate for a mod- 
ified theory of gravity that can be used in the understanding of different astrophysical phenomena, 
usually described by dark matter. The simplest way to do this in the future is by direct compar- 
isons with observed gravitational lenses (particularly Einstein's rings) in clusters of galaxies. We 
are working in this direction and the results obtained will be published elsewhere. 
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